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We argue that coherent production of hard dijets by lin- 
early polarized real photons can provide direct evidence for 
chirality violation in hard processes, the first measurement of 
the magnetic susceptibility of the quark condensate and the 
photon distribution amplitude. It can also serve as a sensitive 
probe of the generalized gluon parton distribution. Numerical 
calculations are presented for HERA kinematics. 

PACS numbers: 12.38.-t, 11.30.Rd, 13.60.-r 



The QCD vacuum is a highly complex state. It has 
a nontrivial particle and energy density characterized by 
quark and gluon condensates and complicated dynami- 
cal properties that characterize its response to external 
probes. In particular, consider quarks and antiquarks in 
the QCD vacuum in a constant (electro)magnetic field 
[|J. In the weak field, the induced magnetization of the 
vacuum is proportional to the applied field, the quark 
density (qq) , the quark electric charge e q and a constant 
X that is called magnetic susceptibility. In relativistic 
notation 



(O\qa a0 q\O) F = e q x(qq)F a0 



(1) 



where F a p denotes the electromagnetic field strength. If 
the magnetic field is varying with a certain frequency, the 
magnetic susceptibility is replaced by the corresponding 
response function which is nothing else but the photon 
distribution function in the infinite momentum frame. 

To be more precise, the wave function of a real photon 
contains both the perturbative chiral-even (CE) contri- 
bution of the quark-antiquark pair with opposite helici- 
ties, and the nonperturbative chiral-odd (CO) contribu- 
tion with quarks having the same helicity and which is 
due to the chiral symmetry breaking. The perturbative 
CE contribution is singular ~ l/|r| at small transverse 
distances r and well known: 



(2) 



iN c e q 
47r 2 r 2 



q+ / du e~ iu ^ (e^ ■ r)(2u - 1) ± ie tk v t e [ k 



where = is the two-dimensional antisymmet 

ric tensor in the transverse plane, e u 
etc. 



The nonperturbative CO 



e d = -(l/3)y / 4ira EM , 
contribution is regular at small transverse separations 
(apart from the logarithms) and can be parametrized by 
the photon distribution amplitude (f>-y(u, /i) Q 

(o\q(o)v a Mx)\j {x) (q)) = 

i 

= ie qX (qq) (e^qp ~ ef\ a ) f due'^ 0>,/x) . (3) 



Here the normalization is chosen as J du(f>-y(u) = 1, and 
u stands for the momentum fraction carried by the quark. 

The distribution amplitude (/> 7 (ii,/i > 1 GeV) is be- 
lieved to be not far from the asymptotic form 



=6tt(l-«). 



(4) 



The magnetic susceptibility was estimated using the vec- 
tor dominance approximation and QCD sum rules 



X (qq) - 40- 70 MeV (at fi = 1 GeV) 



(5) 



However, any direct experimental evidence on both x and 
4>j(u) is absent. The CO contribution in photoproduc- 
tion was discussed only once in j^] for the vector meson 
production at large t. 

In this letter we argue that this structure can be stud- 
ied in experiments for the exclusive hard dijet production 
off nucleons (and nuclei) 



7 + N -> (qq) + N , 



(6) 
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similar to the recent studies of coherent dijet produc- 
tion by incident pions by the E791 collaboration ||. In 
particular, we will show that perturbative (chirality con- 
serving) and nonperturbative (chirality violating) contri- 
butions can be separated by the different dependence on 
the longitudinal momentum of the dijets and on the az- 
imuthal angle. 

The approach of Q is different as compared to to ear- 
lier studies of the dijet photoproduction in that the 
exclusive dijet final state is identified by requiring that 
the jet transverse momenta are compensated to a high 
accuracy within the diffractive cone and making some 
additional cuts. This approach seems to work for the 
case of coherent dijet production from nuclei by incident 



1 



pions, and for photoproduction the corresponding exper- 
imental program is under way at HERA We assume 
that separation of the exclusive qq dijet final state is fea- 
sible. 





FIG. 1. Sample diagrams for the hard dijet photoproduc- 
tion, see text. 

The kinematics of the process and the notation for the 
momenta is shown in Fig. |l|. The Sudakov decomposition 
of the jet momenta with respect to the momenta of the 
incoming particles p\ and P2 reads 



2 2 

Qi - 1\ 

qi = zpi H p 2 + q±,q2 = zpi + ~P2 - q±_ 

zs zs 



(J) 



so that z is the longitudinal momentum fraction and q± 
the transverse momentum of the quark jet. Hereafter we 
use the shorthand notation u = (1 — u) for any longitu- 
dinal momentum fraction u. Note that we consider the 
forward limit, when transverse momenta of the jets com- 
pensate each other. In this kinematics the invariant mass 
of the produced qq pair is equal to M 2 = q 2 ^/ zz, and the 
momentum of the outgoing mic\conp' 2 = P2(l — £)/(1+0j 
where £ = M 2 /(2s ~ M 2 ) ~ M 2 /2s, s = {p x + p 2 ) 2 . 

Since the CE and CO contributions lead to final states 
with different helicity, they do not interfere and the di- 
jet cross section is given by the incoherent sum, for the 
linearly polarized photon 



2tt 



da 



7 — >2 jets 



dtpdq^dtdz 



t=0 



(1 -4zz cos 2 cj>)\JcE\ 2 + 



l aEMogr, r 2 (l 
' 2N c ql 

N 2 q\ 



\Jcof 



(8) 



where cj> is the azimuthal angle between the jet direc- 
tion and the photon polarization (e^ -qx) ~ cos</>, Jce 
and Jce are the CE and CO amplitudes, respectively, 
and the prefactors are introduced for future convenience. 
Note that the CE contribution is ~ 1 jq\ || and the CO 
contribution is suppressed by one extra power of q\ which 
follows from twist counting. The different <fi dependence 
can be traced to the fact that the qq pair is produced 
in a state with orbital angular momentum L z = and 
L z = ±1 for the CO and CE contributions, respectively. 

The CE contribution originates from the region of large 
momenta flowing through the photon vertex and was con- 
sidered previously in [ fL0| p"3| in the high energy limit us- 
ing k± factorization. We believe that the collinear fac- 
torization is more adequate for HERA energies and, in 



difference to the dijet production by incident pions p4[ , 
expect that it is valid for the CE amplitude to all or- 
ders in perturbation theory. Hence the amplitude Jce is 
given by the convolution integral of the coefficient func- 
tion and the generalized parton distribution, cf . jl4| . For 
high energies the gluon contribution is dominant. To 
leading order (LO) in the strong coupling a s = a s (q±) 
the amplitude is given by the sum of Feynman diagrams 
of the type shown in Fig. [l]a with all possible attachments 
of the gluons. The answer reads 



rjLO 

Jce 



(9) 



where H g (y,^) = H g (~y 7 £,) is the generalized gluon dis- 
tribution in the symmetric notation [ jl5| ] , y + £ and y — £ 
are the i-channel gluon momentum fractions with respect 
to (p2 +p' 2 )/2. The full calculation of the NLO contribu- 
tion goes beyond the tasks of this letter. For the reasons 
explained below it is necessary, however, to include the 
leading NLO contribution at large energies (enhanced by 
log£) which corresponds to an additional gluon exchange 
between the i-channel gluons, see Fig. [lp. Including this 
contribution, the imaginary part of the amplitude equals 



a,N r 



dy 



« y + 



« (10) 



where T~t' g (£,£) — dH g (y,£)/dy\ y =£. The real part is 
smaller and can be neglected in the first approximation. 

For high energies alias y — ► 0, H g (y, £) ~ y~ A , where 
in perturbation theory A ~ a s log s jq\ has to be treated 
as a small parameter. Therefore, despite the fact that 
the two terms in ( ^0| ) appear in different orders in the 
collinear expansion, they are of the same order as far 
as the counting of energy logarithms is concerned. This 
feature is specific for real photons and can be traced to 
the fact that the LO amplitude in (J9J) only contains a 
(rather unusual) double-pole, but no single poles (c.f. 
(|j"T|)). Since H g (y,^) ~ G(y) at y ^> £, and as the factor 
oisN c / (ny) is nothing but the low-y limit of the DGLAP 
gluon splitting function, the integral in Eq. ([To]) can be 
identified to logarithmic accuracy with the unintegrated 
gluon distribution /(£, q 2 ) = <9G(£, q 2 )/d\uq 2 . This con- 
tribution corresponds to the one considered in [|l0|-|l2|] in 
the k± factorization approach. The first contribution in 
© is analogous to Eq. (42) in @. 

For the nonperturbative CO contribution the large mo- 
menta are not allowed in the photon vertex and the fac- 
torization formula contains a convolution with the pho- 
ton distribution amplitude. In this case an additional 
hard gluon exchange is mandatory and the diagram in 
Fig. [l]b presents one example of the existing 31 LO con- 
tributions, cf. Fig. 11 in Ref. Jl4| ). The calculation of 
this contribution is similar to the case of pion diffraction 
dijet production considered in much detail in |T3]. Here 
we only present the final result (C F = (N 2 - ljf(2N c )): 
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uu 2N r \u 



(y-£-ie) 



(11) 



Integration over the quark longitudinal momentum frac- 
tion contains a logarithmic divergence at the end points 
U — > 0, 1 which signals that the collinear factorization is 
violated. The divergent contribution is purely imaginary 
and reads 



z 2 + z 2 
2N C 



(12) 



where we have introduced an infrared cutoff u m in = 
Mir/ 1 ?!- I n numerical calculations we use /iir, = 
500 MeV. The origin of the factorization breaking is that 
both final and initial state interactions are present and 
lead to pinching of the integration contour in the so-called 
Glauber region, see |l4|] for a detailed discussion. In the 
present context violation of factorization is probably not 
surprising since the CO contribution is suppressed by a 
power of q\ compared to the leading twist. 

Another important integration region for the quark 
momentum fraction in (O) is £ <C \u — z\ <C 1 which 
produces a logarithmic enhancement in energy: 



rrDLA 



UN c ^ 7 (z) 



dy 



W fl (v,0. (13) 



Hence in the double-logarithmic approximation collinear 
factorization is valid for the CO contribution as well. 

Assuming that the photon distribution amplitude 
</> 7 (z,/i = qj_) is close to the asymptotic form (Q), we 
obtain Jco ~ #(1 — z) for both integration regions, up 
to small corrections. The IR divergence in ( p"2] ) does not 
have, therefore, any significant effect on the jet distribu- 
tion but mainly influences the normalization. 

In the numerical calculation presented below we have 
taken into account the full result in ( |10|) and the imag- 
inary part of the CO contribution in (|11|). Calculations 
are done for HERA kinematics, s = 10000 GeV 2 and the 
value x{qq) = 50 MeV, cf. (||). We use the parametriza- 
tion of the generalized gluon distribution by Freund and 
McDermott JT(| that is based on the MRST2001 leading- 
order forward distribution |]l7f . The transverse momen- 
tum dependence of the cross section integrated over (f>, z 
and t |Q is shown in Fig. [|. 




FIG. 2. The cross section d<T 7 +jv_2jcts+jv/dgi, see text. 
The dashed curve corresponds to the calculation with the 
lower limit in the integral in ( |Io| ) changed to 2£ and a dif- 
ferent IR-cutoff u m = 350 MeV. 

For q± > 4 GeV the cross section is dominated by the 
perturbative CE contribution. The nonperturbative CO 
contribution at q± ~ 4 — 6 GeV is of the order of 



daco/dq\ 
da C E/dq\ 



(7 ± 2 GeV) 2 



x{qq) 

50 MeV 



(14) 



Note a large mass scale which is due to the kinemati- 
cal enhancement of the CO amplitude by a large factor 
4-7r 2 , compare (B) and (||). On top of this, the two con- 
tributions in ( Jl 0[ ) have opposite sign and tend to cancel. 
For large transverse momenta the first term dominates, 
which is natural since it is leading order in the collinear 
expansion. However, progressing towards lower q± one 
effectively goes over to smaller values of Bjorken x ~ 2£ 
so that the second contribution eventually becomes larger 
than the first one and the imaginary part of CE ampli- 
tude changes sign. This explains an abrupt change of 
the slope of the cross section at q± ~ 3 GeV. For smaller 
transverse momenta the dijet cross section is dominated 
by the CO contribution. Because of the cancellations in 
the CE contribution at q± ~ 2.5 — 3.5 GeV the complete 
NLO calculation is required to make the predictions in 
this region fully quantitative. As a crude estimate of the 
uncertainties we show as dashed curve in Fig. ^ the re- 
sults if the lower limit in the integral in ( |To| ) is changed 
to 2£ and with a different IR-cutoff /iir. = 350 MeV. 

The transition between the two different regimes is 
seen very clearly from the dependence of the cross section 
on the dijet longitudinal momentum fraction and the az- 
imuthal angle, compare Fig. |and Fig. @. At q ± = 5 GeV 
the (parton level) z-distribution is almost flat, while the 
cf) distribution is almost purely ~ 1 — cos 2 (f>. In contrast 
to this, at q± = 2 GeV the z-distribution is comparable 
with ~ z 2 (l — z) 2 while the ^-distribution is flat. 

To avoid misunderstanding, we repeat that all calcu- 
lations in this paper are done for the fixed value x(qq) — 
50 MeV and one particular model |l6| of the generalized 
gluon distribution. The related uncertainties are not in- 
cluded. In fact, coherent photoproduction of dijets may 
present the best opportunity to constrain both of them. 
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FIG. 3. The differential cross section 

dcr 7+ jv_2jcts+JV /dq±dz (a) and , da J+N -, 2 jcts+N / 'dq\d(f> (b) 
for jet transverse momentum q± = 5 GeV. Identification of 
the curves is the same as in Fig. |i[ 



the collinear factorization is, strictly speaking, broken. 
However, the sensitivity to the IR cutoff is relatively 
weak and can formally be eliminated by taking into ac- 
count Sudakov-type corrections in the modified collinear 
factorization framework. We think that this technique 
is potentially more accurate and the results can be im- 
proved systematically by the calculations of higher-order 
corrections. In particular, the complete NLO calculation 
of the perturbative CE contribution would be very wel- 
come because of cancellations that are discussed in the 
text. 
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FIG. 4. Same as in Fig. | but for q± = 2 GeV. 

To conclude, we summarize our main points. In this 
letter we argue that studies of exclusive photoproduc- 
tion of dijets with large transverse momenta can yield 
important information on the photon structure at small 
distances. Our main result is that the nonperturbative 
CO contribution is large in the region of intermediate 
q± ~ 2 — 4 GeV and can be clearly separated from 
the perturbative contribution by a different z- and 4>- 
dependence. Observation of the CO contribution would 
be the first clear evidence for the chirality violation in 
hard processes and also provide the first direct measure- 
ment of the magnetic susceptibility of the quark conden- 
sate. On the other hand, the dijet cross section for large 
q± can serve to constrain the generalized gluon distribu- 
tion. 

On the theoretical side, we deviate from previous stud- 
ies of the dijet production by consistently applying the 
collinear factorization in terms of generalized parton dis- 
tributions. For the nonperturbative CO contributions 
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